X-ray reflected spectra from accretion disk models. 
III. A complete grid of ionized reflection calculations 
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ABSTRACT 



We present a new and complete library of synthetic spectra for modeling 
the component of emission that is reflected from an illuminated accretion disk. 
The spectra were computed using an updated version of our code xillver that 
incorporates new routines and a richer atomic data base. We offer in the form of 
a table model an extensive grid of reflection models that cover a wide range of 
parameters. Each individual model is characterized by the photon index F of the 
illuminating radiation, the ionization parameter £ at the surface of the disk (i.e., 
the ratio of the X-ray flux to the gas density), and the iron abundance A-p e relative 
to the solar value. The ranges of the parameters covered are: 1.2 < T < 3.4, 
1 < £ < 10 4 , and 0.5 < A Fe < 10. These ranges capture the physical conditions 
typically inferred from observations of active galactic nuclei, and also stellar- 
mass black holes in the hard state. This library is intended for use when the 
thermal disk flux is faint compared to the incident power-law flux. The models 
are expected to provide an accurate description of the Fe K emission line, which 
is the crucial spectral feature used to measure black hole spin. A total of 720 
reflection spectra are provided in a single FITS filqj suitable for the analysis 
of X-ray observations via the atable model in xspec. Detailed comparisons 
with previous reflection models illustrate the improvements incorporated in this 
version of XILLVER. 



1 http : //hea- www. cf a. harvard. edu/~ javier/xillver/ 
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1. Introduction 

A large number of observations in the past few decades have established the general 
properties of the X-ray spectra from accreting luminous sources. These are systems in which 
a compact object, either a neutron star or a black hole, affects the surrounding material 
with its vast gravitational forces, which leads to accretion of large quantities of gas and 
consequently to the release of large amounts of radiation, mostly in the X-ray band. The 
compact object can be a supermassive black hole as those present in the active galactic 
nuclei (AGN) of many galaxies, or of stellar-mass, such as is the case of many galactic black 
holes (GBHs). 

In AGN, the X-ray continuum is usually characterized by a power-law spectrum 
with photon index T, which can typically be observed to lie in the 1.8 > T > 2.2 range, 
although in some extreme cases it can be either closer to 1 or larger than 3, extending to 
high-energies with an exponential cutoff around 100 — 300 keV. At energies below 0.1 keV, 
the power-law often meets a soft excess that mimics a blackbody radiation. In GBHs, 
the X-ray spectrum in most cases is dominated by a thermal blackbody-like component 
that peaks at high-energies (0.2 < kT < 2 keV), but can also show a high-energy tail 
component of emission. In general, this spectral energy distribution is explained in terms of 
a geometrically thin, optically-thick accretion disk around the compact object. The energy 
dissipation within the disk would be responsible for the quasi-blackbody emission observed. 
The power-law continuum is believed to originate through the Compton up- s cattering o : 



the thermal photons by the electrons in a hot corona or a jet (lHaardt 
1 



1993 



Dove et al 



199 



Dauser et al. 



20131). 



The presence of this dense (rin ^ 10 12 cm -3 ), warm (T ~ 10 5 — 10 7 K), and 
optically-thick (tt > 1) medium is also supported by the detection of atomic features from 
several ions. These and other features constitute an important component of the X-ray 
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spectrum observed from accreting sources, resulting from the reprocessing of radiation by 
the material in the disk. This component is commonly referred to as reflection, in the sense 
that it is the result of radiation that is returned from the accretion disk by fluorescence or 
electron scattering. The current paradigm is that the original power-law radiation irradiates 
the surface of the accretion disk. The X-ray photons then interact with the material 
producing diverse atomic features. These can be produced both via absorption (mostly 
in form of edges), and emission (in form of fluorescence lines and radiative recombination 
continua, RRC). Therefore, the reflection component provides direct information about 
structure, temperature, ionization stage, and composition of the gas in the accretion disk. 

The presence of the Fe K-shell fluorescence emission and the absorption K-edge 
observed in the 6 — 8 keV energy range are recognized as strong evidence for reflection. 
X-ray photons that are photoelectrically absorbed have enough energy to remove a Is 
electron from its K-shell, leaving it in a quasi-bound state above the continuum (autoionizing 
state). The K-hole is then filled by an electron, and the energy difference can be released by 
emitting a second electron (Auger process), or by the emission of a K-shell photon. These 
transitions are of the parity-changing type np — Is, where n = 2 and n = 3 correspond 
to the Ka and K/3 transitions, respectively. Higher n transitions are allowed but with a 
much lower probability. The fluorescence yield, i.e. the probability of emitting a photon 
over an Auger electron, is proportional to the nuclear charge Z to the forth-power (i.e., 
oc Z A ), making the Fe K-shell emission particularly strong. This has bee n shown to be 



true in a large number of observations fi 


1995; 


Winter et al. 


2009; 


Ne et al. 


2010) 



rom X-ray accreting sources (e.g. 



Gottwald et al. 



20101 ) . The Fe K emission line has also proven to be 



crucial in the determina tion of one of the fundamental qu a ntities that describe black ho 



its angular momentum (ILaor 



1991 



Dabrowski et al. 



19971 : 



Brenneman fe Reynolds 



20061 ). 



cs, 



If the reflection occurs near the black hole, line photons will suffer Doppler effects, light 
bending, and gravitational redshift, which produces a skewed line profile with a red wing 
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that can extend to very low en e rgies, particularly in the case of high spin (e.g. . iFabian et al. 



2000 



2008 



Fab ian fe Vaughan 



Steiner et al. 



2011 



2003 



Reynolds fe Nowak 



Reynolds et al 



2012 



2003 



Dovciak et al. 



Dauser et al 



2004 



Miller et al 



20121 ). Therefore, the proper 



modeling of the reflected disk component is of vital importance for determining one of the 
two fundamental parameters that define a black hole. 

The first theoretical studies of X-ray reflection assumed that irradiation on the 
surface of the accretion disk was weak enough so the gas remai ns neutral, but yet wou ld 



reprocess the radiation producing observable spectral features (IGuilbert fe Rees 



1988 1. 



Gr een's functions to des c ribe t h e scattering of photons by cold electrons were first derived 



by 



Lightman fe Rybickil (1198 



observations discussed in 



Lig htman et al. 



Lightman fe White! ( 



1988 



(119811) . and their implications for AGN 



). The latter approach is use d for the 



calculation of cold reflection in the model pexrav ( 
none o f these calculations included line production. 



Magdzi arz fe Zdziarski 



George fe Fabian 



1995 



( jl99lh and 



However, 



Matt ct al. 



( Il99ll ) included the X-ray fluorescence line emission in their Monte Carlo calculations, 



providing line stre ngth, angular distribution, and equivalent widths for the Fe K line. 



Zycki et al. 



(1994) carried out similar calculations including photoionization equilibrium. 



yet neglecting the intrinsic emission inside the gas. 

Much more detailed calculations of the radiative trans 
medium were carried out by 



Ross et al. 



( 11978( 1 and 



er of X-rays in an optically-thick 



Rossi ( I1979I ). Their code solves the 



transfer of the continuum photons using the Fokker-Planck diffusion equation, including a 
modified Kompaneets operator to properly treat the Compton scattering, while the transfer 
of lines is calculated using the escape probabilities appr oximation. This code has been 



updated over the years leading to the reflionx model (IRoss fe Fabian 



1993 



20051 ). which 



has bee n widely used t o mod el the reflected component in the spectra of accreting X-ray 



sources. 



Dumont et al 



(2003) used the TITAN code to examine the accuracy of the escape 
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probability methods versus the "exact" solution of the radiative transfer by imp l emen ting; 
accelerated lambda iterations. This code has been extended by 



Rozahska et al. 



fl2002h to 



treat the cases of Compton-thick media. All these calculations assume constant density in 
the material. It has been argued that a plane-parallel slab under hydrostatic equilibrium 



could represent the surface of an accretion disk more accurately (iNayakshin et al 



20001) 



and that its refle cted spectrum is in fact different from the one predictec 



models (see also 



Rozanska fc Czernv 



1996 



Navakshin fc Kallman 



2001 



)y constant densit y 



Ballantvne et al 



2001 



Dumont et al. 



2002 



Ross &: Fabian 



2007 



Rozahska &: Madej 



20081 ) 



Besides the techniques used to treat the transfer of photons through the media, the 
codes mentioned above may also differ in terms of the atomic data used, which in most cases 
offers a limited perspective of the physics governing the atomic processes for the absorption, 
excitation, and emission p rocesses. These limitations have been over come by our new 



reflection model xillver (IGarcia &: Kallman 



2010 



Garcia et al. 



20111 ). xillver calculates 



the reflected spectrum emerging from the surface of an X-ray illuminated accretion 
disk by simultaneously solving the equations of radiative transfer, energy balance, and 
ionization equilibrium in a Compton-thick, plane-parallel mediu m. The transfe r is solved 



using the Feautrier method under a lambda iteration pro cedure (Mihalas 



1978 



XILLVER 



200lh to 



implements the well-known photoionization code xstar (IKallman fc Bautistal 
solve the ionization structure of the atmosphere, therefore making use of the most updated, 
accurate, and complete atomic database for the X-ray line emission, photoabsorption, and 
photoionization for all the astrophysically relevant ions. 

In this paper we present a complete library of reflection spectra using an updated 
version of our code xillver. This grid of models covers a wide range of input parameters, 
relevant to model the spectrum from accreting X-ray sources. Each model is defined by 
the photon index T, the ionization parameter £ (given by the ratio of the X-ray flux over 
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the gas density), and the abundance of Fe with respect to its solar value. The models are 



provided i n a single tab 



as xspec (jArnaud 



suitable 



;o be used in the commo nly used fitting packages such 



19961 ). and ISIS (IHouck fe Denicola 



20001 ) . We show a detailed analysis 



of our models along the space of parameters and a careful comparison with other similar 
models. The improvements achieved are stressed, and their implications on the analysis of 
X-ray spectra is discussed. 

This paper is organized as follows. In Section [2] we describe the basic aspects of 
the theory used in our calculations, paying particular attention to the changes and 
improvements implemented in the code. An analysis of the new reflection models for the 
different parameters is presented in Section |3j We show the effect of the photon index, the 
ionization parameter, and the iron abundance on the ionization structure and ultimately on 
the emergent spectra. A thorough comparison with other reflection models is also provided. 
The main conclusions and future prospects are summarized in Section HI 



The Reflection Model 



In order to calculate the reflected spectra from X-ray illuminated accretion disks we 
m ade use of our reflection c ode xillver. The details of the calculations are fully described 



in 



Garcia fe Kallmanl (120101 ) . thus we shall review only the main aspects, in particular those 
where some changes have been applied. One important modification takes place in the 
solution of the radiation transfer equation. This equation describes the interaction of the 



radiation field with the gas in the illuminated atmosphere, and it is now expressed as 



Or 2 



u(ji,E,t)-S(E,t) 



(1) 



2- 



http : / /hea-www . cf a . harvard . edu/~ j avier/xillver/ 
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where u(fi, E, r) is the energy density of the radiation field for a given cosine of the incident 
angle with respect to the normal /i, energy E, and position in the slab, which now is 
specified in terms of the total optical depth 



dr = x{E)dz, 



(2) 



where 



X (E) = a kn (E) + a a {E) (3) 

is the total opacity of the gas that includes both scattering and absorption. Here, a>kn{E) 
and a a (E) are the scattering and absorption coefficients, respectively. The former is given 
by the product of the gas density times the Klein-Nishina cross section o~k n {E) for electron 
scattering (which includes the relativistic correction). The latter is defined by the product 
of the gas density times the absorption cross section o~ a (E) due to bound-bound, bound- free, 
and free-free processes. The source function S(E,r) in Equation ([rj, is given by the ratio 
of the total emissivity to the total opacity. It can written as 



S(E,r) = ^lj c (E,r) + mT) 



(4) 



X (E,t) ^ ' ' x (E,ry 
where j(E, r) is the continuum plus lines emissivity, and J C (E, r) is the Comptonized mean 
intensity of the radiation field resulting from the convolution 



J c (E,t) = [ dE'J{E',T)a kn {E')P{E',E). 

0-kn{E) J 



(5) 



Here, J(E, r) = J u(fi, E, r)dfi is the unscattered mean intensity of the radiation field. The 
quantity P(E', E) is the probability of a photon with energy E to be Compton scattered to 
an energy E', which we approximate by assuming a Gaussian profile centered at 

AkT E 



E C = E 1 + 



¥ 



i r 



m P c z 



(6) 



mass et al.lll978[ IRoss fc Fabianlll993l ). with energy dispersion 

1/2 

2k I 2 I h. \~ 
o = E + - 



m P cr 



(7) 
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where k is the Boltzmann's constant, T is the gas temperature, m e is t 
the el ectron, and c is the speed of light. In Equation fl5]) we now follow 



re rest-mass of 



Nayakshin et al. 



( 120001 ) and properly take into account the relativistic (Klein-Nishina) correction to the cross 



section for electron scatt ering, rather than using the classical (Thomson) approximation 



( iGarcfa &: Kallman 



20101). 



Referring to Section 2 in 



Garcia &: Kallman! (120101 ) . the radiative transfer equation 



written in terms of the Thomson optical depth tt has an extra linear term when compared 
with the usual form shown here (Equation [TJ. We found that the radiation transfer equation 
written in terms of the total optical depth behaves better when there are large gradients in 
the opacity, despite the fact that the integration along the path is performed in a different 
grid for each energy bin (i.e., the total optical depth is a function of the photon energy). 
We use a logarithmically spaced grid fixed in tt with 200 spatial bins in the 10 -4 < tt < 10 
range. At each energy, the total optical depth is calculated by integrating Equation (j2J). 



The boundary conditions are only slightly different from those in 



Garcia fc Kallman 



( 120101 ). At the top of the slab (r = 0), we specify the radiation field incident at a given 
angle /i by 



du(r, fi, E) 



u(0, /i, E) 



IF, 



6(n - fio), 



JO ^0 

where F x is the net flux of the illuminating radiation integrated in the whole energy band. 
At the inner boundary (r = T max ), we specify the outgoing radiation field to be equal to a 
blackbody with the expected temperature for the disk: 

du(r, /i, E) 



dr 



+ u(T maX} fi } E) = B(T C 



disk I 



(9) 



Shakura fc Sunyaev 



where B{T) is the Planck function, and Tdisk can be defined using the 
( 119731 ) formulae. Nevertheless, for the models presented in this paper we have chosen to 
neglect any illumination from below the atmosphere, i.e., B(T) = 0. This is convenient 
when comparing our models with previous reflection calculations such as reflionx 
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( IRoss fc Fabianll2005f ). Also, for the parameters typical in AGN, the temperature of the 
disk is sufficiently low that the intrinsic black body emission is weak compared to the 
power-law incident at the surface. Since these models are calculated under the assumption 



of co nstant density, we use the common definition of the ionization parameter (ITarter et al. 



19691 ) to characterize each case, namely, 



£ = — 10 

where F x is the net integrated flux in the 1-1000 Ry energy range, and n e = 1.2nn, and 
n e and nn are the the electron and hydrogen number densities, respectively. The solution 
of the system is found by forward elimination and back substitution. With a solution for 
the radiation field u((jl,E,t), a new J c (E,t) and thus a new S(E,r) can be calculated. 
These are then used to update the solution for u([i, E, r). The procedure is repeated until 
these quantities stop changing within a small fraction. A full transfer solution must be 
achieved iteratively in order to self-consistently treat the scattering process. In general, this 
procedure requires ~ r^ a:r iterations for convergence. Since our calculations are carried out 
up to 10 Thomson depths, we perform 100 iterations. This ensures convergence for energies 
where scattering is dominant, for which tt ~ r. For energies where the photoelectric opacity 
is large, the total optical depth can be much larger than 10. Nevertheless, if absorption 
dominates over scattering, the first term in the right-hand side of Equation (j3J) is reduced. 
In the limit a a ^> «fc n , the source function essentially becomes independent of J C (E, r), and 
the solution converges very rapidly. 

The structure of the gas is determined by solving the ionization balance equations 
for a given gas density and for a particular solution of th e radiation field. At each point 



200l|) to 



within the slab, we use the photoionization code xstar (IKallman fc Bautistal 
calculate level populations, temperature, the total opacity x(£',r), and the emissivity 
j(E,r), assuming that all the physical processes are in steady-state and imposing radiative 
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equilibrium. We have updated our code xillver to work with the current version of xstar 
(version 2.2. lbn) f[ and the most recent atomic database. Although the physics is the same, 
this new version includes improved routines that considerably reduce the computing time, 
allowing the exploration of a wider parameter space for a reasonable allocation of resources. 



2.1. Atomic Data 



As mentioned before, xillver implements the xstar routines for the calculation of 



the ionization structure of the gas, while making use of xstar's atomic d ata. 



the xstar atomic database is described in detail in 



Bautista fc Kallman 



been construct ed using inform ation fr om many sources such as CHIANTI ( Landi &: Phillips 



20061 ). ADAS flSummers 



and the IRON project (IHummer et al. 



2004 ), NIST 



Ralchenko et al. 



he core of 



(2001). It has 



2008|), TOPbase flCunto et al. 



19931 ) 



19931 ). Over the last several years, we have dedicated 



a significant effort towards the investigation of the K-vacancy states in many ions. This 
has allowed us to systematically improve the modeling of the K lines and edges relevant for 
high-quality astronomical X-ray spectra. 

Such atomic data sets include energy levels, wavelengths, Einstein ^-coefficients, 
radiative and Auger widths computed for a large number of ions using three di fferent 



atom ic-structure theoretical a pproaches, inc 



981 ), the autostructure ( 



( jBerrington et al. 



1987 



Badncll 



Seaton 



1986 



uding the relativistic Hartree-Fock ( ICowan 



19971 ). and the multi-configuration Dirac-Fock 



19871 ) methods. Extensive calculations of photoabsorption 
and photoionization cross sections have been also performed using the Breit-Pauli R-Matrix 
method i ncluding the effects of radiative and Auger damping by mean s of an optical 



potential ( iGorczyca fc Robicheaux 



1999 



Gorczyca fc McLaughlin 



2000J). 



3- 



http : //heasarc . gsf c . nasa . gov/xstar/xstar . html 
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In particular, calculations of the atomic data required for the spectral m odeling of the 



jC-she ll pho 



toabsorption of oxy gen and nitrogen ions have been carried out in 



Garcia et al 



2005|) and 



Garcia et al. 



( 120091 ). Computations have also been carri ed out for the atomi c 



structu re of the isonuclear s equences of Ne, Mg, Si, S, Ar, and C a (IPalmeri et al. 



nickel (IPalmeri et al. 



2008bl ). and aluminum (IPalmeri et al. 



2008ah . 



201ll ). More recently, we have 



focused our attention on the iron peak and light odd-Z elements, i.e., F, Na, P, CI, K, Sc, 
Ti, V, Cr, Mn, Co, Cu and Zn producing the atomic parameters for more than 3 million 



fine-structure K lines (IPalmeri et al. 



20121 ). Photoabsorption calculations of the cross 



sections across the K-ed ge of Ne, Mg, Si, S, Ar , and Ca have been performed for ions with 



less than 11 electrons by 



Witthoeft et al. 



Witthoeft et al. 



(120091 ). and for ions with more than 10 electron by 



(l2011bl ). In these calculations it was shown that damping processes affect 
the resonances converging to the K thresholds causing them to display symmetric profiles of 
constant width that smear the otherwise sharp edge at the photoioniz a tion th resholds. The 



Witthoeft et al. 



( l2011aT ). These new 



Li-like to Ca-like ion stages of nickel are discussed in 
data sets are continually being incorporated into the xstar data base in order to generate 
improved opacities in the K-edge regions of the ions considered. 



Th e main differences in the current version of the atomic data with respect to the one 



used in 



Garcia fe Kallmanl ( 120101 ) are as follows: The inclusion of the atomic data for the 



n 



2—7-3 inner-shell transitions of Fe vi-xvi ions, the so-called Fe M-shell unresolved 



transition array (IGu et al 



20061 ). The inclu sion of the K-shell at o mic data for the 



S, Ar, Ca, Al. and Ni isonuclear sequences 



Palmer i et al. 



2011 



and dielectronic recombination rates from 



(IPalmeri et al. 



2008a 



Mg^Si 



Witthoeft et al. 



2011b 



Witthoeft et al.l 12011a ). The implem entat ion of the updated r adiative 



Badnelll (I2006h and 



Badnell et al. 



(120031 ). Also, 



we have improved the resolution of the high-energy extrapolation implemented in the Fe 
photoionization cross sections for energies well above the K-edge. Finally, by looking at 
the third- row Fe ions (Fe I-Vlll), we found that for the densities used in our simulations 
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(n e ~ 10 15 cm -3 ), many metastable states are populated. Previously, we only used K-shell 
photoionization data for the ground state, therefore missing the contribution from these 
metastable states. This resulted in a weak Fe K emission line in models with low ionization 
parameter. This problem has been corrected in the new version of the data base by 
including the photoionization cross sections for all these states. 



2.2. Illumination spectrum 



Another important change in the xillver code is the definition o 



Garcia fe Kallman 



;he radiation field 



fl2010h . the incident 



that illuminates the surface of the accretion disk. In 
spectrum was assumed to be a power-law with a given photon index T for all energies in 
the 10 -1 — 2 x 10 5 eV range. This kind of illumination is appropriate for spectra with 
r ~ 2. However, for steeper spectra (T > 2) this choice of illumination creates unphysical 
conditions, since the incident power-law will tend to have very many photons concentrated 
in the low energy part of the spectrum. Thus, we have now adopted a power-law that breaks 
at 100 eV and decays with an exponential tail for lower energies. We have also extended the 
energy to 1 MeV, and moved the exponential high-energy cutoff from 200 keV to 300 keV. 
Figure [T] shows the impact of this modification for the resulting reflected spectrum in a 
case where T = 3. The solid lines show the reflected spectra predicted by xillver using 
the old (black) and the new (red) definitions of the illumination. The incident power-law 
is shown with dashed lines for each case. Both calculations have a similar ionization 
parameter (log £ ~ 2) but the integrated flux is different given the differences in the shape 
of the power-law. Thus, the model with the new definition (broken power-law) needed 
to be scaled down such that both incident spectra will have the same flux near 100 eV 
and energies above, region in which the photon index is T = 3. The huge impact of the 
illuminating spectrum on the reflected component is clear, especially at energies below 
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5 keV. This difference becomes more dramatic as the photon index increases. However, 
for harder spectra (T < 2), the break in the incident power-law has little impact on the 
reflected spectrum, especially in the high-energy band. 

The extension of the spectrum to higher energies has also an important influence on 
the reflection calculations. This is particularly important in models on the other extreme 
of the photon index parameter space. In the case of a very hard illumination spectrum 
(r < 2), models with similar ionization parameters will differ in the number of photons 
that are concentrated in the high energy part of the spectrum. Figure [2] shows the resulting 
xillver models for an incident power-law with T = 1.4, log £ = 2.8, and solar abundances. 
For the black and red curves, the energy range extends upward to 200 keV and 1 MeV, 
respectively. In panel (a), the dashed lines show the incident power-law spectra, while solid 
lines are the reflected spectra for each case. The emission lines in the model with the high 
energy coverage (red) are considerably weaker than those in the original model (black). 
This is most evident for the Fe K complex around 6 — 7 keV, and for the O VII La near 
650 eV. The reason for this difference is that in models with the high-energy extension, 
the region of the atmosphere near the surface is more strongly heated by the high-energy 
photons, and thus its temperature is higher than in the original model. This can be seen 
in panel (b) of Figure [2} where the temperature profiles in the vertical direction of the disk 
are plotted as function of the Thomson optical depth. Clearly, the two models converge at 
large depths (tt ~few), but the one in which the illumination extends to higher energies 
results in a considerably hotter atmosphere near the surface. Although we expect the Fe 
features to be produced well inside the disk, line photons that escape through a hotter gas 
will tend to suffer a larger energy shift from Compton scattering, with similar probabilities 
for up- or down-scattering. This makes the line profile broader and more symmetric. 
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3. Results 

As in the previous version of xillver, all the models shown here are calculated 
with an energy spectral resolution comparable with the best current observations 
(1Z = E/AE ~ 350), which requires 5000 bin points in the energy range considered here 
(10 _1 — 10 6 eV). Each model covers a large column density of gas (tt = 10), using 200 
spatial zones and 10 bins that describe the angular dependence of the radiation fields. 
The first step in the spatial grid has been reduced to tt ~ 10~ 4 to ensure an accurate 
representation of the illuminated layers in the slab. The hydrogen number density is held 
constant at n-^ = 10 15 cm -3 for all the models presented in this paper. 



3.1. The space of parameters 



A new library of reflection spectra has been produced using the new version of the code 
xillver. This set of models covers a wide range of parameter values that are relevant for 
fitting the spectra of accreting sources. Each model is characterized by the photon index 
r of the incident power-law radiatio n, the ionization parame ter £, and the iron abundance 



Ap c with respect to the solar value (IGrevesse fc Sauval 



19981 ). In order to have a library of 



models suitable for both AGN and GBHs, we have produced models covering photon indices 
in the T = 1.2 — 3.4 range, in steps of 0.2. Observations typically show that the reflection 
signatures from accretion disks in AGN are produced in a material at a lower ionization 



stage than those observed from GBHs (IGarcia et al. 



201ll ). consistent with the picture 



that the accretion disk in the latter are much hotter than in the former. Consequently, we 
vary the ionization parameter over a wide range to cover both classes of sources. We have 
produced models with £ = 1, 2, 5, 10, 20, 50, 10 4 , 2 x 10 4 , 5 x 10 4 erg cm s -1 . Finally, we 
also allowed the iron abundance to be treated as a free parameter, given the importance of 
the Fe K emission profile in most accreting sources. For simplicity, the abundance of all 
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the other elements conside red, namely H, He, C, N, 0, Ne, Mg, Si, S, Ar, and Ca, is kept 



fixed to the solar values of 



Grevesse fe Sauvall (119981 ). Thus, models with A Fe = 0.5, 1, 5 and 



10 were calculated, taking into account both sub- and super-solar Fe abundances. Here, 
A-pe = 1 corresponds to an iron abundance of 2.5 x 10~ 5 with respect to hydrogen. These 
choices of parameter values resulted in a library of 720 synthetic reflection spectra that can 
be used for the modeling of the reprocessed component in X-ray observations. The present 
set of models is provided as a single file □ in FITS (Flexible Image Transport System) 
format, which can be loaded into the fitting package xspecJ^I via the atable model. 



3.2. The effect of varying the ionization parameter £ 

Figure [3] shows a sub-group of the resulting reflected spectra for three values of the 
photon index T and a range of ionization parameters. The Fe abundance is set to the 
solar value (A Fe = 1) for all the models shown here. Panels (a), (b), and (c) correspond to 
T = 1.4, 2 and 2.6, respectively. Each panel shows the reflected spectra calculated for a 
different value of the ionization parameter. From bottom to top, each curve corresponds 
to £ = 1, 2, 5, 10, 20, 50, 100, 200, 500, 1000, 2000, 5000 erg cm s" 1 , respectively. Because 
the gas density is held fixed at n e = 1.2 x 10 15 cm -3 for all these calculations, increasing 
the ionization parameter is equivalent to increase the flux F x of the illuminating source, 
where F x = £n e /47r. The spectra are plotted in units of EF E (equivalent to vF v if plotted 
in frequency), so that a power- law of T = 2 is shown as a horizontal line. Note that no 
rescaling or shift is applied; instead, each curve is color-coded according to its corresponding 
value of log £, to improve clarity. This Figure provides a general overview of the effect that 

4 http : / /hea-www . cf a . harvard . edu/~ j avier/xillver/ 



5- 



http : //heasarc . gsf c . nasa . gov/ xanadu/xspec/ 
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the ionization parameter has on the spectrum reflected from an optically thick plane-parallel 
slab. The original power-law shape of the illuminating continuum (for the lowest value of 
£), shown in black dashed lines, suffers drastic modifications due to both absorption and 
emission. The spectra are shown in the entire energy range included in the calculations 
(0.1 — 10 6 eV). In the low-energy part of the spectrum (0.1 — 10 eV), the continuum is 
dominated by bremsstrahlung emissivity, despite the fact that this is where the illumination 
flux is decreased due to the cut-off imposed at 100 eV. Emission lines and absorption edges 
from H, He, and C are clearly visible for most of the models with low and intermediate 
ionization. In the 10 — 10 4 eV energy range, the photoelectric opacity dominates over 
the electron scattering opacity. Therefore, this region is where most of the absorption 
occurs, yielding large departures from the original power-law continuum. Bremsstrahlung 
emissivity decreases rapidly for energies above ~ 100 eV; however, many emission lines 
from all the ions included in these calculations remain visible over the entire energy range. 
At higher energies (> 10 4 eV), electron scattering is the dominant source of opacity since 
the cross section for photoelectric absorption decays as ~ E~ 3 , while the Klein-Nishina 
cross section for electron scattering remains fairly constant. Thus, in this spectral region 
Compton scattering is the only relevant process and the reflected spectrum depends on the 
shape of the original illuminating field. 

By looking at the overall shape of the reflected continuum shown in Figure [31 it is 
possible to distinguish between two main regimes: (1) The high-ionization case, where the 
resulting spectra mostly shows very narrow emission features, while the continuum still 
resembles the original shape of the illuminating power-law; and (2) the low-ionization case, 
where the emerging spectra are a combination of a very rich and complex set of emission-line 
profiles superimposed on a strongly absorbed and modified continuum, which significantly 
departs from the original power-law. The specific value of the ionization parameter that 
separates these two regimes depends, to a certain degree, on the photon index V. In fact, 
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for the T = 1.4 case shown in panel (a) this transition is quite obvious, as a drastic change 
in the reflected spectra can be seen between models with log £ = 2.3 and 2.7. These changes 
are the result of large differences in the ionization balance solutions, as can be seen in 
Figure HJ which shows the corresponding temperature profiles for each one of the models 
shown in Figure [HI color-coded in the same way. Because larger ionization implies a larger 
illumination flux, models with high £ values are systematically hotter than those with low 
£ values. The illuminating radiation, incident at the surface of the slab (at tt = 10 -4 for 
our proposes), penetrates the first layers heating the gas. The amount of heating not only 
depends on the illuminating flux, but also on the particular shape of the radiation field; 
e.g., note that the models for small values of T are consistently hotter than those with large 
values. We shall discuss these effects further in the next Section. As photons are scattered 
or absorbed and re-emitted at different energies, the shape of the continuum is modified 
reducing the net heating in the gas. At some point H and He recombine which increases the 
cooling very rapidly. The gas then suffers a rather sudden transition to a lower temperature, 
as can be seen for most of these models. This transition occurs at deeper regions in the slab 
for larger ionization parameters and for softer input spectra (large T). Regardless of the 
ionization parameter or the shape of the ionizing radiation, in all cases the field thermalizes 
and the gas eventually reaches the same lower temperature around (0.5 — 1) x 10 5 K. This 
is the limit in which the illuminating radiation no longer contributes to the ionization of 
the material and its temperature is simply set by the density of the gas. 

To illustrate in more detail the effects of the ionization parameter on the reflected 
spectra, we will concentrate on models with one particular value of the photon index. 
Further, we will examine the low- and high-ionization models separately. Thus, Figure |5] 
shows some of the reflected spectra resulting from models with T = 2 in the 10 2 — 5 x 10 5 eV 
energy range, which is the spectral band typically covered in X-ray observations. Panel (a) 
shows 4 models with low-ionization, specifically £ = 1,5,20, and 100, multiplied by factors 



-20 - 



of 1, 10 2 , 10 4 , and 10 6 , respectively. Their corresponding temperature profiles are curves 
1, 3, 5 and 7, from left to right, in panel (b) of Figure |H The dashed line represents 
the illuminating power-law spectrum for the model with the lowest ionization parameter 
(£ = 1). As mentioned before, these models display a strong decrease of the continuum flux 
due to photoelectric absorption, in particular for energies below ~ 20 keV. Nevertheless, a 
very rich and complex set of fluorescent emission lines due to K-, L-, and M-shell transitions 
from many ions is also present, superimposed on a highly absorbed continuum. The general 
spectral shape of the lowest ionization models (£ = 1 — 20) is very similar, where only small 
changes are seen. The temperature in these models is relatively low (T < 10 6 K), and the 
gas settles to a low temperature regime at a small optical depth (tt < 0.1), meaning that 
most of the slab remains neutral. The emission due to the Ka and K/3 transitions in Fe is 
distinctive at ~ 6.4 and ~ 7.1 keV, respectively. The case for £ = 100 is where the changes 
in the reflected spectra become more evident, in both the continuum absorption and in 
the emission features. The higher temperature in this model allows for the excitation of 
elements such as Mg, Al, Si, S, Ar and Ca, which produces more fluorescence lines in the 
1 — 10 keV region. The Fe Ka emission becomes slightly broader and the K/3 is less intense. 
Some radiative recombination continua (RRC) can also be seen near ~ 0.5 — 1 keV. 

Panel (b) of Figure shows the reflected spectra for the next five values of the 
ionization parameter considered here (£ = 200, 500, 1 x 10 3 , 2 x 10 3 , and 5 x 10 3 ). Each curve 
has been rescaled by a constant factor to improve clarity. These factors are 1, 10, 10 2 , 10 3 , 
and 10 5 from low- to high-ionization, respectively. As before, the dashed line shows the 
illuminating power-law spectrum corresponding to the model with the lower ionization 
(£ = 200). These are the models considered as high-ionization, since the continuum is not 
highly affected by the photoelectric opacity. This can be seen by comparing, for example, 
the reflected flux at 10 2 eV in the spectrum for £ = 200 and in the one for £ = 1, shown 
in panel (a). As mentioned before, increasing the ionization raises the temperature of the 
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illuminated region of the slab. This region also extends to deeper zones as the radiation 
is able to ionize the gas at larger optical depths, as can be seen in the corresponding 
temperature profiles (see panel b of Figure H]). The increase in the temperature has two 
main effects. On the one hand, it affects the ionization of the gas. The ions from low-Z 
elements are completely stripped from all their electrons, while the ions of the heavier 
elements are partially ionized. This changes the emission lines produced inside the gas, 
and thus the emerging spectra lacks emission from low-Z ions, progressively showing more 
emission from highly ionized O, Ne, Ar, Ca and Fe, as the illumination increases. On the 
other hand, the high temperature affects the way line photons are scattered by electrons, 
given the dependence of Equations ([6]) and (J7|) on the kinetic energy of the electrons. When 
4kT ~ E, the probability for a photon to either gain or lose energy after each scattering 
becomes comparable, which effectively produces a broadened and more symmetric line 
profile. This effect is particularly evident in the Ka emission lines from O and Fe, observed 
at ~ 0.65 and ~ 6.9 keV, respectively. We shall return to a more detailed discussion of the 
spectral features in Section 13.51 

3.3. The effect of varying the photon index T 

It is clear from the discussion in the previous section that both the net flux and the 
spectral shape of the ionizing radiation incident on the surface of the slab have a great 
impact on the ionization balance of the gas, and thus on the reflected spectrum that will 
emerge at the surface. For this particular library of models we have adopted a power-law 
shape for the illumination spectrum, thus the general shape is controlled by changing the 
value of the photon index T. Figure E] shows a sub-group of the reflected spectra calculated 
for various conditions. Panels (a), (b), (c) and (d) show the resulting models for a given 
ionization parameter, i.e., £ = 10, 10 2 , 10 3 , and 10 4 , respectively. Each panel contains the 



-22 - 



spectra of models for all the values of the photon index F considered in our calculations 
(T = 1.2 — 3.4), color-coded accordingly. Since the goal of this Figure is to show the general 
trends in the spectra introduced by changing T, all the curves are plotted with the same 
normalization (thus, no rescaling was applied). It is quite obvious how the changes in the 
ionizing continuum affect the resulting continuum of the emergent radiation, as expected. 
However, the ionization structure of the gas seems to follow different trends on T for 
different values of £. 

By comparing the two extreme values of £ (panels a and d), a completely opposite 
behavior is seen. In the low ionization case (£ = 10), models with low T values are the 
ones with the most absorption, suggesting a much colder gas than those for high T values. 
This is consistent with the temperature profiles obtained for these models, which are shown 
in panel (a) of Figure [71 Therefore, in the low-ionization regime, the softer the input 
spectrum (large T value), the hotter and more ionized the gas becomes. Looking at the 
high- ionization case in panel (d) (£ = 10 4 ), the opposite occurs: the models with harder 
illumination spectra (lowest T) are completely ionized, owing to the large gas temperature. 
Thus, in the high-ionization regime, the harder the input spectra, the hotter and more 
ionized the gas becomes. The intermediate regimes plotted in panels (b) and (c) of Figures [6] 
and [7J show the transition between these two regimes. For £ = 10 2 the gas temperature is 
low for the harder spectrum (r = 1.2), and increases with the photon index up to T ~ 2.4, 
where it starts decreasing again. For the case with £ = 10 3 (panel c), the transition occurs 
at T ~ 1.8. 

The reason for this change of behavior is related to the processes that contribute 
to the heating and cooling of the gas. In general terms, there are two main competing 
mechanisms: photoionization heating plus recombination cooling, and Compton heating 
and cooling due to electron scattering. The low ionization models are those for which the 
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illumination is relatively low. In this case, photoionization is the dominant process that 
heats the illuminated layers of the slab. Because the heating rate due to photoionization is 
essentially given by the radiation field flux times the photoelectric opacity of the gas, and 
since the latter is dominant at energies below ~ 10 keV, a very hard input spectrum will 
produce much less heating than a soft spectrum, owing to the lack of low energy photons. 



The contrary is true in the high-ionization regime where the dominant proce ss is Compton 

Garcia fc Kallmanl ( 2010 ). in this limit 



heating and cooling. As discussed in Section 3.1 of 
the gas temperature approaches an asymptotic value, the Compton temperature, given by 

< E > 



T, 



where 



c 



< E >-- 



4k 



J F(E)EdE 



(11) 



(12) 



J F(E)dE 

is the mean photon energy, which is a quantity that only depends on the spectral shape. In 
Figure M we show the resulting Compton temperature as a function of the photon index 
T in the range of our calculations, which agrees very well with the temperature of the hot 
layer in the high-ionization models shown in panel (d) of Figure 0J Physically this makes 
sense, since for the input spectra with low F values most of the photons are concentrated in 
the high energy part of the spectrum, where Compton scattering becomes very important. 
Therefore, in the high-ionization regime, hard spectra are more efficient in heating the 
illuminated layers of the slab. 



3 Note, however, that xstar (and consequently xillver) employs a 



ment of the Compton heating and cooling from 
Equation (fTTj) . 



Guilbert et al. 



fill relativistic treat- 



( 119821 ) . rather than using 
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3.4. The effect of varying the Fe abundance 

The elemental abundances considered in a photoionization calculation can potentially 
affect the ionization balance, temperature structure, and ultimately the observable spectral 
features in the reprocessed radiation. The total amount of a particular element changes the 
continuum opacity, which in turn affects the photoionization heating rate. At the same 
time, the abundance of a particular element influences the strength of the emission and 
absorption features due to bound-bound and bound-free transitions. Given the relevance 
of the Fe emission in the analysis of the X-ray spectra from accreting sources, we have 
carried out calculations in which the Fe abundance, normalized to its solar value, is 
varied between sub-solar, solar, and super-solar values. All the other elements considered 
in these calculations are set to their solar values. Figure [9] shows a comparison of the 
reflection calculations for different values of the iron abundance A-p e . Left panels show the 
temperature profiles, while right panels show the corresponding reflected spectrum in the 
f — 10 5 eV energy range. In all the panels, each curve corresponds to one particular value 
of Af p = 0.5, 1, 5 and 10 , where A Fe = 1 corresponds to 2.5 x 1CT 5 of Fe with respect to 



H ( iGrevesse et al. 



19961 ). In each one of the right panels, the plotted spectra have been 
rescaled for clarity. The scaling factors are, from bottom to top, 10~ 2 , 1, 10 2 , and 10 4 . 
Top, medium, and bottom panels correspond to ionization parameters £ = 10, 10 2 and 10 3 , 
respectively. The photon index is set to T = 2 in all these models. The general tendency is 
the same in all these simulations. The increase of the Fe abundance induces more heating in 
the illuminated layers due to the increase in the opacity (and thus a larger photoionization 
rate), which raises the gas temperature. However, because continuum absorption is also 
increased, the radiation field thermalizes at a smaller depth for the high abundance models, 
as can be seen from the temperature curves. Note that the increase in temperature is more 
subtle in the high- ionization case (top panel), since in this regime the Compton heating and 
cooling is the dominant process that controls the gas temperature. 
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The effects of the Fe abundance are evident in the reflected spectra as well, both in 
the emission and the absorption features. The 10 2 — 10 4 eV energy range clearly shows a 
substantial reduction of the flux due to the increase in the continuum opacity; meanwhile, 
the Fe K edge near 8 keV grows deeper as Ap e becomes larger. At the same time, all the 
Fe emission features are affected as well, which is mostly evident in the Fe K emission 
complex in the 6 — 8 keV region. The strength of the whole emission profile increases when 
the abundance is high, as expected. In the high ionization models (top panel), there is 
a distinctive RRC profile right before the Fe K edge, which becomes substantially more 
prominent for Ap e = 10. Also, there are enhanced emission lines near 100 eV and just 
above 1 keV, which correspond to M- and L-shell transitions in iron, respectively. However, 
these emission profiles are noticeable only for the £ = 100 model (middle panel), in 
particular the L-shell. The reason for this is that these transitions occur in a rather narrow 
range of ionization stages. If the gas is very neutral, the photoelectric absorption reduces 
considerably the number of ionizing photons at those energies, reducing the number of 
excitations from the L-shell. On the contrary, if the gas is hot and ionized, the fraction of Fe 
ions with L-shell electrons is very low (i.e., Fe i-xvi), as most of them are stripped already. 
This is of great relevance, since this feature can be used to constrain both the ionization 
of the gas as well as the iron abundance. One particular example is the Seyfert 1 galaxy 
1H 0707-495, which X-ray spectrum shows evidence of a very intense emission in both the 
Fe K- and L-shell regions. Fits using reflection models re quire a high Fe abundance at an 



ioniz ation parameter consistent with the present analysis (IFabian et al. 



20121 ). 



2012 



Dauser et al 



Although somewhat extreme, 1H 0707-495 is not the only case where Fe is found to be 
over-abundant based on predictions from reflection models. In fact, at least for AGN, it is 



common 



spectra (IFabian 



y the case th at super-solar iron abundance is required to fit the observed X-ray 



20061 ). A possible explanation for the apparent extreme Fe abundance 
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in some AGN has been proposed by 



Reynolds et al 



(120 12l ). on the grounds of radiative 



levitation in the accretion disk. If the radiative force exerted on a Fe ion exceeds the vertical 
gravity, this could cause iron to diffuse towards the photosphere of the disk, enhancing its 
abundance. But, in general, there is no particular reason why the other elements should be 
considered at their solar values, except the simplicity inherent in such an approximation. 
Given the capability of xillver to treat any particular choice of elemental abundances, it 
can be used to produce smaller set of models custom made for any specific situation. This 
could be of great use in peculiar systems such as ultra-compact X-ray binaries, where a 
prominent O VII Ka emission line observed in the X-ray spectra is thought to originate from 



reflection in an accretion disk over-abundant in oxygen (Madej et al. 



2010 



Madej fe Jonker 



201ll ). Also, emission lines from H-like S, Ar and Ca ions detected in some lo w-mass X-ra y 



binaries suggest reprocessing material with 



2009 



di Salvo et al 



2009 



Egron et al 



20JJ). 



compositions different from solar (ID'Ai et al. 



3.5. Spectral features 

Figure [3] shows the great complexity of the reprocessed spectrum emerging from 
an illuminated, optically-thick slab. Meanwhile, the models presented here constitute 
a high-resolution representation of only one of the components observed in the X-ray 
spectrum of accreting sources, namely the reflected component. In reality, one observes a 
composite spectrum that includes the original power-law (presumably coronal) component 
plus possibly a thermal blackbody-like component. Additionally, if the reflection occurs 
within a few gravitational radii from the compact object, Doppler and gravitational redshifts 
will smear the spectral profiles. Finally, absorption due to intervening gas such as warm 
absorbers and outflows can also be present. We shall discuss some of the most prominent 
and representative of these features that one expects to observe given the capabilities of 
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current detectors. 



3.5.1. The Fe K-shell emission 

Undoubtedly, the emission complex near 6 — 7 keV, which is due to transitions from 
the inner K-shell of Fe ions, is the most prominent atomic feature in the X-ray spectrum 
of accreting sources. It is this feature that provides the clearest evidence for reflection of 
high-energy photons in the relatively cold, optically-thick material of an accretion disk. The 
ubiq uity of this feature has been established observationally f or a large number of sources 



(e.g. 



Gottwald et al. 



1995 



Winter et al. 



2009 



Ng et al. 



2010 



Fukazawa et al. 



201 lh for the 



following two reasons. First, the fluorescence yield, i.e. the probability of emission of a 
photon rather than an Auger electron, is proportional to Z 4 , where Z is the nuclear charge. 
Second, the Fe K-shell lines are emitted in a clean region of the X-ray spectrum (6 — 8 keV), 
where few other ions emit or absorb radiation. Furthermore, in this energy range galactic 
absorption is negligible and most detectors operate quite effectively. 

Figure [10] shows a compilation of all the radiative transitions from Fe ions in the 



Kallman et al 



6 — 1 keV energy range that are included in our databasqj(cf., Figure 3 in 
2004 ). There are a total of 2735 lines, all of which correspond to K-shell transitions. 
The open circles show the line energy plotted against the ionization stage of each ion, as 
determined by the relation Z — N e + 1, where Z is the nuclear charge and N e is the number 
of electrons in the ion. Thus, an ionization stage of 2 corresponds to Fe II (single ionized), 3 
to Fe ill (double ionized), and so forth. Filled circles show the most intense transitions, i.e., 
those with large transition probabilities (here we have chosen A r > 10 13 s _1 ). Ionization 
increases upward: The H- and He-like ions are near the top and the neutrals near the 



7 http : //heasarc . gsf c . nasa . gov/uadb 
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bottom. The big group of points to the left corresponds to the Ka transitions (n = 1 — > 2, 
with n being the principal quantum number), and the smaller group to the right (for which 
most energies are above 7 keV) corresponds to the K/3 transitions (n — 1 — Y 3). Note that 
what is commonly referred to as the neutral Fe K line at 6.4 keV is in fact a combination of 
several transitions that span the energy range 6.39 — 6.43 keV; this feature can be produced 
by many ions ranging from Fe I up to Fe xvil. For Fe xvill and more ionized ions, the 
transition energy spread is larger and the average energy moves monotonically towards 
higher energies as the ionization stage increases. For the He- and H-like ions (Fe xxv-xxvi), 
there are only a few lines with energies around 6.9 keV. Notice also that the K/3 transitions 
are only produced up to Fe xvil, since for higher ionization stages the n = 3 shell is empty. 
The line energy for the K/3 transitions varies between ~ 7.04 keV for Fe II up to ~ 7.19 keV 
for Fe XVIL 

It is these line energies and intensities that we use to analyze the emission in the Fe 
K region of the reflected spectra. Figure [11] shows, in the 5 — 8 keV band, the reflected 
spectra for T = 2 and Ap e = 1 for different values of the ionization parameter in the range 
£ = 5 — 5 x 10 3 . Ionization increases downward in the figure. Starting at the lowest 
level of ionization (£ = 5 — 50), the first few models show similar groups of Fe lines, with 
two distinctive emission lines centered at ~ 6.4 keV and ~ 7.1 keV, which correspond 
respectively to the Ka and K/3 lines. The line energies indicate that the emission is 
dominated by Fe ions with ionization stages lower than Fe xvil. The smaller feature below 
6.4 keV, which is produced by Compton down-scattering of line photons, is usually referred 
to as the Compton shoulder. For models with £ > 200, the K/3 line is no longer visible, 
which suggests that the emission is dominated by ions at higher ionization stages. In 
particular, models with £ = 200 and 500 show a rich complex of emission lines at energies 
between 6.4 — 6.7 keV. At £ = 10 3 , the emission is centered at ~ 6.7 keV and ~ 6.9 keV, 
which implies that the gas is highly ionized and most of the Fe ions are hydrogenic or 
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He-like. In this state the temperature is high, and the overall line emi ssion profile is 



therm ally broadened into a symmetric profile via Compton scattering. 



Kallman et al 



( 120041 ) present a more detailed but similar analysis of the Fe K emission complex based 
on xstar simulations, which includes atomic data for Fe; however, they do not consider 
energy redistribution due to Compton scattering. 



3.5.2. Lower-Z elements 

In addition to the important Fe-K emission complex in the reflected spectrum, our 
code includes many emission lines of astrophysically relevant ions at lower energies, which 
are due to such elements as C, N, O, Ne, Mg, Si, S, Ar, and Ca. The intensities of these 
lines tend to increase as the illuminating spectrum softens (i.e., as T increases) and as the 
flux correspondingly increases at lower energies where photoelectric absorption dominates. 
Figures [T21 d31 and IT41 show the reflected spectra for T = 3 and A Fe = 1 and for all the values 
of the ionization parameter considered in our library. The spectra cover the 0.3 — 4.5 keV 
energy range, which is the band that contains most of the inner-shell transitions from low-Z 
ions. At the top of each figure is indicated the energy of the strongest Ka emission line for 
each ion in the isonuclear sequence of each of the elements contained in our atomic database 
(except for the neutral and single- ionized cases). While these Ka lines are not at all the 
only lines considered in our models, they nevertheless give an indication of the energy range 
one expects most of the features of a given element to appear. 

Figure [121 shows spectra for the lowest ionization models, namely £ = 1, 2, 5, 10, and 20 
(with ionization increasing from bottom to top). At the lowest energies (0.3 — 1 keV), the 
lines are indistinguishably blended together. At higher energies, emission lines from low 
ionization states of Mg, Si and S are evident, and relatively weak lines due to Ar and Ca 
Ka are also present. The radiative recombination continuum (RRC) due to H-like Ne at 
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~ 1.4 keV becomes visible in models with £ > 5. Such features are produced by recombining 
electrons with energies that exceed the ion binding energy. The excess energy is radiated 
as a photon. The energy of a typical RRC photon is E = Ejp + kT e , where Ejp is the 
ionization potential of the ion and T e is the electron temperature. 

For higher ionization parameters, the reflected spectra are dominated by RRC features, 
as shown in Figure Q3J which displays spectra ranging from £ = 50 at the bottom to £ = 10 3 
at the top. Also visible are higher-excitation emission lines of Mg, Si, S and Ar. The H-like 
oxygen RRC at ~ 0.87 keV is quite strong for £ = 50 (bottom curve), and it weakens as 
£ increases. Strong emission features due to Ka transitions of Mg, Si, and S appear at 



1.4, 1.8, and 2.45 keV, resp ectively. The 



Ka line and the Si XIII RRC (IGarcia et al. 



after feature corresponds to a blend of the S XV 



201l|). 



Figure [TH shows the last five models of the series with the ionization parameter 
increasing from £ = 2 x 10 3 (bottom) to £ = 5 x 10 4 (top). For these large values of £, fewer 
emission lines are present. However, they are quite intense because the soft illuminating 
spectrum is rich in low-energy photons. The strongest emission lines are produced by H-like 
ions of O, Mg, Si, and S with energies of ~ 0.654, 1.47, 2.01, and 2.63 keV, respectively. The 
RRC from H-like Si and S are clearly seen at ~ 2.67 and 3.5 keV, respectively. Generally, 
Ka fluorescence from the low ionization states of O, Mg, Si, S, Ar, and Ca ions is most 
important at low column depths because of the higher K shell opacity of these elements 
relative to Fe. Conversely, Fe Ka is more important at relatively high columns, i.e., Tt ~ 1. 



3.6. Comparison with previous models 



Considering the several reflection models currently available, reflionx (IRoss fe Fabian 



20051 ) is the one that is most similar to the models presented in this paper. It has been 
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widely used by the scientific community in analyzing spectral data for many observations of 
various black-hole and neutron-star sources. We therefore benchmark our results by making 
a detailed comparison of xillver and reflionx, while describing the advantages of our 
model. 

In reflionx, the radiation field is separated into the direct component of the 
illuminating radiation, which is unaffected by either scattering or absorption (assumed 
to be oc e _TV ), and the diffuse component, which results from Compton scattering and 
emission within the slab. The Compton processes are treated by sol ving the Fokker-Planck 



diffusion equation, which includes a modified Kompaneets operator (IRoss et al. 



1978 



Ross 



19791 ). Incident photons that are Compton scattered in the slab contribute to the emissivity 
of the diffuse field, and it is assumed that their distribution in energy is described by a 
Gaussian with central energy and energy dispersion given by Equations ([6]) and (JTj). Due 
to the deficiency of the Fokker-Planck equation in handling steep gradients in the energy 
spectrum, resonance lines a re treated using the es cape probabilities technique. In addition 



to fully ionized species, the IRoss fe Fabianl (120051 ) calculations include: C ill-iv, N m-vii 



O ill- vill, Ne ill-x, Mg iii-xii, Si IV- xiv, S iv-xvi, and Fe iv-xxvi. However, none 
of the neutral or near-neutral ions are include d in their models. The cr oss sections for 
photoionization are calculated from the fits of I Verner fe Yakovle v ( 1995 1) ; in the case of 



Kaastra fe Mewel f ll993[ ) 



Fe, transition probabilities and Auger decay rates are taken from 
The items we have just highlighted are among the most important differences between the 
atomic data employed in the two models (see Section I2TTT) . It is also important to mention 
that the Fe Ka lines treated in reflionx are the recombination lines of Fe xxvi and 
Fe xxv, and the fluorescence lines of Fe vi-xvi, while the Ka fluorescence of Fe xvil-xxil 
is assumed to be suppressed by autoionization. All K/3 and higher n resonances are also 
neglected. 
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The solar abundances of the elements in reflionx are taken from Morrison fe McCammon 



1983), while in xstar, and thus in xillver, we adopt the values of 



Grevesse fe Sauval 



19981 ). Table [T] shows the values reported from both sources for each element, and the last 



column shows the ratio of the two values. Apart from N and Ar, all the abundances used in 
xillver are lower than those used in reflionx. In particular, O is lower by ~ 10% and 
Fe by ~ 30%, while, most notably, Ne is ~ 80% lower. The choice of abundance model is 
relatively unimportant, as illustrated in Figure [T5l which shows reflected spectra computed 
using xi llver for £ = 10 (le f t pan el) and £ = 10 3 (right panel). Results obtained using the 



model of 
of 



Grevesse fe Sauvall ( 



higher 



Morrison fe McCammon 



1998) are shown in black, and those obtained using the model 



Morrison fe McCammon! (I1983T) in red. The differences in each panel are small. For the 



119831 ) abundances, the continuum is slightly depressed and 
some of the strongest lines are somewhat weaker. These small differences are unlikely to be 
important in analyzing real data. 

Although the solar abundances assumed for both reflionx and xillver are fixed and 
different, both models allow the Fe abundance to be varied. This allows us to compensate 
for the abundance differences in making direct comparisons of the two models. One such 
set of comparisons is shown in Figure [161 where the emergent spectra for different values of 
both the ionization parameter and the photon index are plotted for xillver (in black), and 
reflionx (in red). Top, middle and bottom panels are for T = 1.4, 2, and 2.6, respectively. 
Each panel shows 3 pairs of curves, corresponding to £ = 10, 10 2 , and 10 3 , from bottom to 
top, respectively. The spectra are from the FITS tables that are accessible via the atable 
model in XSPEC. To ensure that the energy resolution is the same for all the spectra, in 
all cases we used a logarithmic dummy response of 1000 energy points over the energy 
range 0.1 — 1000 keV. Furthermore, appropriate grid point were chosen in the parameter 
space to avoid problems that might be introduced by the interpolation procedure. The 
normalization for both models is the same because we use consistent definitions for the 
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ionization parameter and the power-law spectrum. Thus, all the parameters were set to 
be identical for the two models, with the exception of the iron abundance: It was set to 
Ap e = 1.32 in xillver in order to compensate for the difference in the solar values assumed 
in each model, as discussed above. 

In general, the two models are in better agreement for low ionization and soft spectra, 
as in the models for T = 2.6 shown in the bottom panel of Figure [TBI However, even in these 
cases there are important differences: For the bottom spectra with £ = 10, the xillver 
model is more absorbed at low energies than the reflionx model. While the energies and 
intensities of the strongest lines are in good agreement, some weaker features at ~ 2 keV 
and ~ 3 keV are only present in the xillver spectrum, probably due to the differences 
in the atomic data. These particular features are much stronger for the harder spectra 
shown in the middle and top panels. The differences in the atomic data also affect the Fe 
K emission profile near 6 — 8 keV. For example, Fe Ka emission at ~ 6.4 keV is somewhat 
more intense in our model and K/3 is conspicuously absent in the reflionx spectra. 
Concerning the Fe continuum, one might expect significant differences because the models 
use very different photoionization cross sections. However, the models agree very well on 
the depth of the Fe K edge near 7.2 keV, especially at low ionization. As mentioned above, 
the discrepancies are most pronounced for high ionization and for the hard spectra shown 
in the top panel (r = 1.4). The difference for the pair of curves for £ = 10 3 is remarkable: 
The xillver continuum below ~ 1 keV is at least an order of magnitude fainter than the 
reflionx continuum; the extremely intense O Lya emission line at ~ 0.65 keV in our 
spectrum is scarcely present in the reflionx spectrum; and the Fe K line in the latter 
spectrum is much broader than in the former. 

These large discrepancies for high ionization and hard spectra are unexpected and 
there is no ready explanation. The lack of emission lines from \ow-Z elements suggests 
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that (for the same input parameters) the slab is hotter and more ionized for reflionx 
than for xillver. We test this idea by comparing the performance of the two models for 
different values of the ionization parameter: In Figure [T71 we compare the two T = 1.4 
models where for reflionx we have doubled the normalization and reduced the ionization 
parameter to £ = 500 (half the value for xillver), in order to match the fluxes at the higher 
energies. The shape and intensity of the Fe K line is now in better agreement. However, 
the reflionx continuum is still considerably higher at low energies; as a consequence, the 
reflected spectrum is even softer than the incident power-law spectrum. 

In order to test whether our code is reliably modeling this low- energy continuum, we 



20121 ). 



compared our results with those computed using the apec model ({Foster et al. 
For this comparison, with xillver we computed the reflection spectrum from a thin slab 
(tt = 10~ 2 ) for solar abundances at very low ionization with the temperature held fixed 
at 10 6 K. These parameters closely match the case of a collisionally-ionized gas whose 
spectrum is dominated by bremsstrahlung emission. Figure [18] shows the comparison. 
Setting aside the differences in the emission lines (a result of the differences in the atomic 
data sets), the spectra are in good agreement, especially the level of the continuum. This 
comparison gives us confidence in xillver's implementation of the free-free emissivity. 

An important quantity to co mpare these models is the equivalent width (EW) of the 



Fe K emission complex. Following 



Garcia et al. 



(1201 if ), we use the well-known formula, 



EW mizl^l dE , (13) 

J 5.5 keV fc{&) 

where F(E) is the flux of the reflected spectrum and F C (E) is the flux in the continuum. 
The 5.5 — 7.2 keV range of integration covers all Fe emission features, including the K/3 lines. 
The continuum is approximated by a straight line that passes through the endpoints of the 
energy band. Clearly, this straight-line continuum and this particular choice of integration 
limits are somewhat arbitrary and by no means constitutes an accurate determination of 
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a physical quantity, but it suffices to compare the two models. Figure [19] shows plots of 
the EWs vs. the ionization parameter for both xillver and reflionx, with F increasing 
from top to bottom. For both models, the EWs tend to decrease as the ionization increases. 
At low ionization, the EWs are similar, although the values for xillver are consistently 
higher, which is expected because the reflionx model lacks the Fe K/3 lines. Meanwhile, 
at high ionization (£ > 10 3 ), where the Fe K emission is dominated by the H- and He-like 
ions, the agreement in the EWs is also good, but values for the reflionx model are in this 
case somewhat greater for large F (bottom panel). 

Significant discrepancies appear when one compares the two models at higher levels 
of ionization, 10 2 < £ < 10 3 : The EWs for the reflionx models decrease drastically, 
especially for the softer spectra (bottom panel). This behavior may occur because reflionx 
does not include the Fe Ka lines for most of the second-row ions (Fe xvii-xx h) under 



the a ssumption that these lines are suppressed by Auger resonant destruction (IRoss et al 



19961 ). xillver, on the other hand, predicts strong Ka emission at the energies at 
which many of these ions are expected to emit (e.g., see Figures fTUl and fTTi) . xstar, and 
consequently xillver, automatically takes into account the effects of Auger decay by 
accurately calculating the branching ratios for Auger versus fluorescence emission; however, 
it does not include resonant scattering of the lines. This omission seems reasonable because 
it is plausible that velocity gradients or turbulence in an accretion disk will remove photons 
from the line core, thereby reducing resonant scattering. Nevertheless, because we disregard 
this mechanism we may be over-predic ting the line intensity for some models. As one 



additional complication, iLiedahll (l2005f ) has shown that Auger resonance destruction is a 
selective process that highly attenuates only a limited subset of Fe Ka lines. We conclude 
that further analysis is required in order to accurately quantify the emission from L-shell 
Fe ions. 
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We have also co mpared our results to t 



lose o btained using the reflection models 



pexrav and pexriv (IMagdziarz fc Zdziarski 



19951 ). which compute the reflected spectrum 
from a completely neutral disk and an ionized disk, respectively. These mod els are an 



19881 ) that models 



implementation of a semi-analytic Green's function (ILightman fe White! 
Compton reflection of X-rays by cold electrons while including bound-free continuum 
opacity. Fluorescence emission line are however not included. The left panel of Figure [2U] 
compares the reflected spectrum predicted by xillver with that predicted by pexrav 
for an incident power-law spectrum with T = 2 and an exponential cutoff at 300 keV. 
Because pexrav models a completely neutral gas, we use the lowest ionization parameter 
considered in xillver, namely £ = 1. The two models agree fairly well in the level of the 
reflected continuum, in particular near the Fe K line region. The edge energy and EW of 
the feature at ~ 7.2 keV are similar. Our model shows more absorption than pexrav in 
the low-energy part of the spectrum. However, the xillver model can not be considered 
completely neutral; there is some heating near the surface (see the leftmost curve in the 
middle panel of Figure Hj) . 

The right panel in Figure [20] compares a reflected spectrum from an ionized disk 
(£ = 10 3 ) computed using xillver with one computed using pexriv. At the outset, we 
note that pexriv, unlike both xillver and reflionx, does not solve for the ionization 
balance of the slab. Instead, pexriv assumes an isothermal medium with a maximum 
temperature of T = 10 6 K. Both the ionization parameter and the gas temperature are 
allowed to vary. This approach results in the large discrepancies that are evident in both the 
level of the continuum and the edge features. The level of the continuum is similar only for 
energies below 1 keV. For higher energies, the pexriv spectrum is highly absorbed and its 
Fe K edge is much deeper. This comparison demonstrates for the case of an ionized reflector 
the importance of accurately modeling the ionization balance, pexriv's assumption of an 
isothermal gas is a poor one for large values of the ionization parameter. Another limitation 
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of pexriv that was pointed out by iFabian fe Rosa (120101 ) is the lack of blurring effects of 
Compton scattering, which smears the K edges. 

We now discuss in turn the differences at high energies (> 20 keV) between our results 
and those delivered by reflionx and pexrav/pexriv. At these energies, photoelectric 
absorption is negligible, and the only source of opacity is Compton scattering. Furthermore, 
the gas temperature has only a very small effect on the Comptonized spectrum. Therefore, 
any differences in the model spectra can be attributed to how each model treats scattering 
of photons by cold electrons, xillver implements a Gaussian convolution for all energies 
and at all depths (Section [2j). reflionx uses this approximation for the lines only, while 
using a modified Kompaneets operator to treat the diffuse radiation. 

We tested the accuracy of our approximation by comparing with results obtained using 



a Monte Ca rlo (MC) co de, w 



the work of 



Matt 



rich we treat as the go 



(J2002|) and 



Magdziarz fe Zdziarski 



d standard. This code, which extends 



(119951 ). simula tes Compton scatteri ng 



and radiative transfer through neutral gas for a semi-infinite slab (jEikmann et al. 



20121 ). 



Included in this code are photoabsorption for all elements with Z < 30 and fluorescent line 
emission using the fluorescent yields from iKaastra fe Mewd (119931 ) . Compton scattering is 



modeled using the proper relativistic scattering dynamics and the differential Klein-Nishina 
cross section. 

Figure [2T] compares reflected spectra computed using xillver and the MC code for 
r = 2, £ = 1 and solar abundances. The result for reflionx is also shown for these same 
parameters. Overall, there is a good agreement between the three models. At energies 
< 10 keV, xillver and the MC simulation predict very nearly the same level of the 
continuum and the same spectrum of intense emission lines. For example, the ~ 6.4 keV 
Fe Ka line has a very similar intensity and shape, although the MC code predicts a 
~ 7.2 keV K/3 line that is somewhat more intense. The main differences occur above the 
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~ 7.2 keV Fe K edge, in the continuum feature referred to as the Compton hump. The MC 
simulated spectrum has a somewhat shallower Fe K edge than the other two spectra, and 
its ~ 20 — 60 keV is somewhat fainter. Above ~ 100 keV, the reflionx spectrum is the 
better match to the MC spectrum. 

As a bottom line, the differences between the MC simulation and either xillver or 
reflionx are not larger than 30%, which is reasonable given the simplistic approximation 
used in those models. However, this level of performance may not be adequate for the 



analysis of such current and fu ture X-ray missions as NuSTAR ( Harrison et al. 



eROSITA flPredehl et al 



20101) and 



20111 ). Furthermore, the Compton-hump spectrum can provide 



important constraints on disk inclination. Therefore, our goal is to improve this aspect of 
the code in future versions of xillver. 



4. Conclusions 



In this paper we have presented a new and complete library of synthetic spectra for 
modeling radiation that is reprocessed in an accretion disk and emitted as a reflected 
spectrum, which is generated in response to illumination by an incident power-law 
spectrum. This pr e sent v ersio n of our code xillv er is an update of those presented in 



Garcia fe Kallmanl (120 lOf ) and 



Garcia et al 



( 120111 ). We have made several improvements to 
both the routines and the atomic data and have produced a large grid of reflection models 
covering a wide range of parameters. Each model is characterized by the photon index T of 
the illuminating radiation, the ionization parameter £ at the surface of the disk, and the Fe 
abundance with respect to its solar value. A total of 720 reflected spectra are provided in 
a single FITS filqj suitable for the analysis of X-rays observations via the atable model in 



http : / /hea-www . cf a . harvard . edu/~ j avier/ xillver/ 
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xspec. In order to represent the physical conditions typically observed in most accreting 
sources, this library covers the following range of parameters: 1.2 < F < 3.4, 1 < £ < 10 4 , 
and 0.5 < A Fc < 10. 

The spectrum that illuminates the surface of the accretion disk is assumed to be a 
power-law in the 10 _1 — 10 6 eV energy range, with a sharp low-energy cutoff at 100 eV 
and an exponential high-energy cutoff at 300 keV. The low-energy cutoff is important for 
spectra with steep power laws (T > 2); it prevents the spectrum from being unphysically 
over-populated with low-energy photons. The power-law is extended up to 1 MeV because 
in the low-r case this portion of the spectrum is rich in high-energy photons, which 
importantly effect both the gas temperature and the reflected spectrum. 

The intensity of the illuminating spectrum, which is specified by the ionization 
parameter £, significantly impacts the ionization structure of the gas. In all models, the 
gas temperature in the vertical direction has a similar profile. Near the surface, where 
the illumination is intense, there is a hot layer (T > 10 6 K) where Compton heating and 
cooling dominate. At larger depths, a warm/cold region exists where photoelectric opacity 
and recombination dominate. The depth at which the transition between these two regimes 
occurs increases as £ increase and as the incident spectrum softens (large T). Above the 
transition layer, the radiation field thermalizes and the gas temperature remains fairly 
constant at around 10 5 K. For low- ionization (£ < 100 — 500), the reflected spectrum 
displays a very rich set of emission lines superimposed on a strongly absorbed continuum. 
For high- ionization (£ > 500 — 1000), the spectrum consists of very narrow emission lines 
from ionized species and a continuum that resembles the incident power-law. 

We have present ed detailed co mparisons with the reflection models reflionx 



(jRoss &: Fabian 



20051 ). and pexrav (IMagdziarz fc Zdziarski 



19951 1 . The reflionx model 



is more similar to ours because of its complexity and the range of physical processes 
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considered. Furthermore, reflionx and xillver both characterize the illuminating 
spectrum as a power law, cover the same energy range, and use the same definition for the 
input parameters. The spectra generated by the two codes are in good agreement. The 
spectra agree best for large T (soft spectra) and low values of £. However, at low energies 
reflionx generates an excess of soft flux compared to xillver. For hard spectra, T = 1.4, 
and high ionization, £ = 10 3 , the difference in the continuum at 1 keV is about an order of 
magnitude. This large discrepancy is not well understood and requires further analysis. 

At very high or very low ionization, xillver and reflionx predict similar EWs for 
the Fe K emission feature (including both Ka and K/3), although the values for xillver are 
generally larger. However, at intermediate levels of ionization, in the 10 2 < £ < 10 3 range, 
the two models disagree strongly: the Fe K EWs computed in this range computed using 
reflionx drop drastically. The effect is presumably because reflionx assumes that the 
Fe K lines from second-row ions are completely suppressed by Auger resonant destruction. 

Comparing XILLVER with the more simplistic models PEXRAV and PEXRIV highlights 
the importance of performing detailed calculations which take into account key physical 
processes. The greatest deficiency of these models is that they do not generate emission 
lines. The neutral pexrav model is in reasonable agreement with xillver concerning 
the level of the continuum and the strength of the edges. However, for moderate or high 
ionization (£ > 20), the pexriv model is in strong disagreement b ecause it assumes an 



isothermal disk, as previously pointed out by 



Fabian fe Rossi (1201 



The library of models presented in this paper is suitable for modeling the reflection 
spectra of accreting sources when the thermal disk component of emission is small compared 
to the incident power-law component. That is, the present version of xillver is suitable for 
analyzing the spectra of AGN, and also GBHs in the hard state, w hen the disk component is 



cool and faint. For GBHs in the thermal or steep power-law states (jRemillard fe McClintock 
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20061 ). the strong thermal component entering the atmosphere from below will significantly 
change the ionization structure of the disk. In a future publication, we will report on an 
extension of xillver that is appropriate for modeling reflection spectra in the presence of 
a strong thermal component. 

Regardless of how accurate xillver (or any such slab reflection model) becomes, it 
is by itself inadequate for modeling the reflection component in the spectra of accreting 
black holes. For example, current black hole spin determinations based on fitting of the Fe 
K line naively assume a single ionization state for the reflecting portion of the disk, which 
can extend over hundreds of gravitational radii. We plan to construct more realistic models 
for particular illumination geometries that allow for an ionization gradient in the radial 
direction. This work will combine t he g eneral relativistic approach implemented in the 



relline code (iDauser et al. 



2010 



20131 1 with the model spectra presented in this paper. 
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Table 1. Elemental Solar abundances 



Element xiLLVER a REFLlONX b Ratio 



H 

He 

C 

N 


Ne 
Mg 

Si 

S 

Ar 
Ca 
Fc 



1.0 
0.1 

3.7 x 10~ 4 
1.1 x 10~ 4 

6.8 x 10~ 4 
2.8 x 10~ 5 
3.5 x 10~ 5 

3.5 x 10~ 5 

1.6 x 10~ 5 
4.5 x 10~ 6 
2.1 x 10~ 6 
2.5 x 10~ 5 



1.0 




1.00 


0.1 




1.00 


x 10~ 


-4 


0.83 


x 10" 


-5 


1.21 


x 10~ 


-1 


0.92 


x 10" 


-1 


0.20 


x 10~ 


-5 


0.88 


x 10- 


-5 


0.94 


x 10- 


-5 


0.84 


x 10- 


-6 


1.18 


x 10- 


-6 


0.94 


x 10- 


-5 


0.76 



Grevesse fc Sauvall (119981 ) 



Morrison fc McCammonl (119831 ) 
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Fig. 1. — Results using different illumination spectra with xillver for a gas with log £ = 
2. The dashed lines are the two different illumination power-law with T = 3. In black, 
the power-law covers the 0.1 — 2 x 10 5 eV energy range, with an exponential cutoff at 
200 keV. In red, the power-law has a sudden cutoff at 0.1 keV, extends to 1 MeV, and an 
exponential cutoff at 300 keV. The black and red solid lines are the resulting reflected spectra 
correspondent to each case. 
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Fig. 2. — Results using different illumination spectra with xillver for a gas with log £ = 2.8. 
In panel (a), the dashed lines are the two different illumination power-law with T = 1.4, 
only differing by their high-energy limit. In black, the power-law covers up to 2 x 10 5 eV 
with an exponential cutoff at 200 keV. In red, the power-law extends to 1 MeV with an 
exponential cutoff at 300 keV. The black and red solid lines are the resulting reflected spectra 
correspondent to each case. Panel (b) shows the respective temperature profiles within the 
slab. 
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Fig. 3. — Reflected spectra for different values of T and £. Panels (a), (b), and (c) correspond 
to T = 1.4, 2.0 and 2.6, respectively. In each panel, solid lines show the reflected spectra 
calculated for a different value of the ionization parameter. From bottom to top, each curve 
corresponds to f = 1,2,5,10,20,50,100,200,500,1000,2000,5000 erg cm s" 1 , color-coded 
accordingly. Note that no rescaling is applied. The illuminating power-law for the lowest 
value of P is shown in black dashed lines. The Fe abundance is set to the solar value ( = 1) 
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Fig. 4. — Temperature profiles corresponding to the models shown in Figure El Panels (a), 
(b), and (c) correspond to T = 1.4, 2.0 and 2.6, respectively. In each panel, the solid lines 
represent the gas temperature along the vertical direction of the slab, corresponding to a 
model with different value of the ionization parameter, color-coded accordingly. The slab is 
illuminated at the surface (tt = 10~ 4 ), and the calculation is carried out up to tt = 10. 



- 54 - 




Energy (eV) Energy (eV) 



Fig. 5. — Results showing the impact of the ionization parameter in the reflected spectra. 
In panel (a), solid lines show the models for £ = 1,5,20, and 100, multiplied by factors of 
1, 10 2 , 10 4 , and 10 6 , respectively, to improve clarity. The dashed-line represents the incident 
power-law for the model with £ = 1. In panel (b), solid lines are the reflected spectra for £ = 
200, 500, 1000, 2000, and 5000, multiplied by factors of 1, 10, 10 2 , 10 3 , and 10 4 , respectively. 
The dashed-line is the incident power-law for the £ = 200 case. In all cases, T = 2 and 
A Fe = 1. 
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Fig. 6. — Results showing the impact of the photon index in the reflected spectra. Pan- 
els (a), (b), (c) and (d) show the resulting models for a given ionization parameter, i.e., 
£ = 10, 10 2 , 10 3 , and 10 4 , respectively. Each panel contains the spectra of models with all 
the values of the photon index T considered in our calculations (T = 1.2 — 3.4), color-coded 
accordingly. No rescaling is applied, and Ap c = 1 in all cases. 
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Fig. 7. — Results showing the impact of the photon index in the temperature profiles. 
Panels (a), (b), (c) and (d) show the resulting models for a given ionization parameter, 
i.e., £ = 10, 10 2 , 10 3 , and 10 4 , respectively. Each panel contains the gas temperature along 
the slab for the models shown in Figure El covering all the values of the photon index V 
considered in our calculations (r = 1.2 — 3.4), color-coded accordingly. Ap c = 1 in all the 
cases. 
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Fig. 8. — Compton temperature as a function of the photon index F for the illuminating 
power-law spectra used in our calculations (see Equation fTTT) . 
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Fig. 9. — Results showing the impact of the Fe abundance. Left panels show the temperature 
profiles, and right panels the corresponding reflected spectrum. Each curve corresponds to 
one particular value of Ap e = 0.5, 1, 5 and 10. In each of the right panels, the plotted spectra 
have been rescaled for clarity. The scaling factors are, from bottom to top, 10 -2 , 1, 10 2 , and 
10 4 . Top, medium, and bottom panels correspond to ionization parameters £ = 10, 10 2 and 
10 3 , respectively. In all these models T = 2. 



-59 - 



+ 

CD 



N 



26 
22 
18 



g> 14 

CCS 

CO 
g 

"cC 

N 

o 



10 
6 
2 





1 1 1 1 I 


1 1 1 i 


1 1 1 i 


■ | i i i i 


i i i i 
® 


i 








o o® 


®® 


msm> 










« 1 ■ t 310 


3 








BID 


' »»» osms 

HO dK 










GO 




© QB 












• 








@> 

O 














<m <s> 














m 














m 


































C3 
O 








m 
<m 




3D) 4 








GDC 
















a» 














mm 














mm 














am 














<m 




.... i 


. . . i 


. . . i 


. i . . . . 


.... 


i 



6.0 6.2 6.4 6.6 6.8 7.0 

Line Energy (eV) 



7.2 



7.4 



Fig. 10. — Emission lines from all the Fe ions in the 6 — 10 keV energy range (open circles). 
In the x-axis is the transition energy, and in the j/-axis the ionization stage of each ion, 
given by Z — N e + 1, where Z is the nuclear charge and N e the number of electrons. Filled 
circles show transitions with the highest probability (A r > 10 13 s _1 ). Data accessible via 



http : / /heasarc. gsf c .nasa. gov / uadb / 
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Fig. 11. — Reflected spectra in the Fe K region (5 — 8 keV). Each panel corresponds to a 
different ionization parameter, as indicated. In all these cases, T = 2 and Ap e = 1. 
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Fig. 12. — Reflected spectra from models with steep illumination (T = 3), and solar abun- 
dances (Ap e = 1). From bottom to top, each curve corresponds to £ = 1,2,5,10, and 20, 
respectively. The strongest Ka emission lines are indicated for each ion considered in the 
calculation (except for neutral and single ionized ions). 
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Fig. 13. — Reflected spectra from models with steep illumination (T = 3), and solar abun- 
dances (Ap c = 1). From bottom to top, each curve corresponds to £ = 50, 100, 200, 500, and 
10 3 , respectively. The strongest Ka emission lines are indicated for each ion considered in 
the calculation (except for neutral and single ionized ions). 
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Fig. 14. — Reflected spectra from models with steep illumination (T = 3), and solar abun- 
dances (A-pe = 1)- From bottom to top, each curve corresponds to £ = 2 x 10 3 , 5 x 10 3 , 10 4 , 2 x 
10 4 and 5 x 10 4 , respectively. The strongest Ka emission lines are indicated for each ion 
considered in the calculation (except for neutral and single ionized ions). 
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Fig. 15. — Ref l ected spectra from xillver us ing the elemental abundanc es of 



Grevesse fc Sauvall ( 119981 ) (black curve), and those by iMorrison fc McCammonl ( 119831 ) (red 
curve). Left panel shows the case for £ = 10, while right panel shows the £ = 10 3 models. 
In all cases T = 2. 
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Fig. 16. — Comparison of the reflected spectra calculated with xillver (black curves), and 
with reflionx (red curves). Top, middle, and bottom panels show the models for T = 1.4, 2, 
and 2.6, respectively. In each panel, pairs of curves correspond to £ = 10, 10 2 , and 10 3 , from 
bottom to top. All spectra are plotted with the same energy resolution. Note that Ap e = 1.32 
in xillver to compensate for the difference in the solar values with those used in reflionx. 
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Fig. 17. — Reflected spectra for T = 1.4 and solar abundances. In black is the xillver 
result for £ = 10 3 . In red is the reflionx spectra for £ = 500 and twice the normalization. 
Although the Fe K emission agrees better, there are large discrepancies in the reflected flux 
at lower energies. 
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Fig. 18. — Comparison of the emergent spectra from a collisional plus bremsstrahlung domi- 
nated gas. In red is the spectrum from the APEC model with T = 10 6 K and solar abundances. 
In black is the xillver prediction for a thin slab (r T = 10 -2 ), at the same temperature, and 
log £ = —5. The blue dashed line shows the incident spectrum. 
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Fig. 19. — Fe K emission line equivalent width versus ionization parameter obtained from the 
spectra calculated by xillver (filled circles), and reflionx (open squares). Top, middle, 
and bottom panels correspond to models with V = 1.4, 2, and 2.6, respectively. 
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Fig. 20. — Comparison of the reflected spectra calculated with xillver (black curves), 
and pexrav/pexriv (red curves). Left panel shows the completely neutral case modeled 
by pexrav, compared with a xillver model with £ = 1. Right panel shows the ionized 
reflection case, thus a pexriv model with £ = 10 3 is compared with xillver at the same 
ionization parameter. The disk temperature in the pexriv model was set to T = 10 6 K (the 
highest possible value). In all cases T = 2, the high-energy cutoff is at 300 keV, and solar 
abundances are used. 
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Fig. 21. — Comparison of the reflected spectra as calculated with xillver (black curve), 
reflionx (blue curve), and the Monte Carlo simulation (red curve), for an illumination 
with r = 2, £ = 1, and solar abundances. 



